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In  recent  years  it  has  become  possible  to  use  density-functional  total  energy 
calculations  to  predict  a  variety  of  properties  of  materials.  Calculations  have 
been  performed  to  predict  the  lattice  constant,  cohesive  energy  and  bulk  moduli 
of  solids1,  ^  ;  high  pressure  phases  of  materials®’  4  ;  geometries,  binding  energies 
and  vibration  frequencies  of  molecules5  ;  surface  energies6  ;  and  the  geometry  of 
monolayer  films7  and  of  surfaces.8,9  Most  total  energy  calculations  of  surfaces 


have  been  confined  to  clean  semi-conductor  surfaces,  the  reason  being  that  the 


pseudo-potential  method  used  in  these  calculations  is  better  suited  to  semi¬ 
conductors  than  to  transition  metals,  either  clean  or  adsorbate  covered.  Very 
recently  there  has  been  some  work  on  clean  transition  metal  surfaces  using  all¬ 
electron  methods.6, 9 


Id  this  letter  we  present  total  energy  calculations  for  an  adsorbate  on  a  tran- 
siton  metal  surface.  The  particular  system  we  have  chosen  is  P(lxl)  H/Ni(100) 
because  a  variety  of  experiments  that  we  can  compare  with:  thermal  desorption10 
,  electron  energy  loss  spectroscopy11  ,  ion-channeling12  and  He  scattering13  have 
been  performed  on  this  system.  We  have  performed  calculations  to  determine  the 
geometry  of  the  adsorbed  H,  its  binding  energy,  its  vibration  frequency  and  the 
barrier  to  surface  diffusion. 

Our  calculations  were  performed  using  the  Linearized  Augmented  Plane 
Wave  (LAPW)  method.14, 15  The  core  electrons  (Ni  ls-3p)  were  treated  fully  rela- 
tivistically  and  the  valence  electrons  scalar  relativistically.  The  core  is  included 
in  the  calculation  of  the  total  energy.  The  potential  was  allowed  to  have  a  gen¬ 
eral  form  everywhere  (no  muffin-tin  approximation),  except  that  non-spherical 
terms  within  the  muffin-tin  spheres  were  evaluated  approximately.  The  non¬ 
spin-polarized  form  of  the  Vosko-Wilk-Nusair  (VWN)16  exchange-correlation 
potential  was  used.  Some  of  the  calculations  were  redone  using  the  Hedin- 
Lundqvist  (HL)17  potential  in  order  to  check  the  sensitivity  of  the  results  to  the 
choice  of  the  exchange-correlation  potential.  Most  of  the  calculations  were  per¬ 
formed  on  7  layer  slabs  consisting  of  5  layers  of  Ni  with  a  layer  of  H  on  each  side. 

Fig.  1  shows  the  adsorption  energy  Eti,  as  a 

function  of  the  vertical  distance  of  the  H  atoms  above  the  surface  Ni  plane  for 
the  center(C),  bridge(B),  and  atop(A)  sites.  Here  E$fNi  and  E^h  are,  respec¬ 
tively,  the  total  energies  (per  unit  cell)  of  the  H  covered,  and  clean  Ni  slabs,  and 
Effem  is  the  energy  of  a  free  spin-polarized  H  atom  (0.957  Ry  or  13.02  eV  for 
VWN,  0.976  Ry  or  13.28  eV  for  HL).  The  equilibrium  position  is  found  to  be  0.6 

o 

a0  (a,  =  0.529  A  )  above  the  center  site  (corresponding  to  a  Ni-H  bond  length 
of  1.79  A).  This  is  probably  a  small  underestimate  («  2%  or  0.07  a,  in  the 
bond  length,  or  0.31  a,  in  the  vertical  distance  z )  since  local  density  functional 


calculations  for  bulk  Ni  gives  a  bond  length  that  is  about  2%  too  small.  Hence, 
we  expect  the  true  position  of  the  H  atoms  to  be  at  about  0.9  at  above  the  sur¬ 
face  Ni  plane.  Also  the  inclusion  of  magnetism  (which  we  have  neglected)  would 
tend  to  increase  interatomic  distances  slightly.  In  fact,  a  recent  spin-polarized 
calculation18  gives  z  =  0.69  a0 .  It  is  difficult  to  determine  the  position  of  H  on 
surfaces  experimentally,19  but  there  is  in  fact  a  reliable  determination  of  the  posi¬ 
tion  of  H  on  Ni(100)  from  transmission  ion-channeling  experiments.12  These 
experiments  show  that  H  adsorbs  0.95  ±  0.2  at  above  the  center  site,  in  excellent 
agreement  with  our  estimate  (0.6  +  0.3  a„ ). 

The  calculated  value  of  the  adsorption  energy  (3.61  eV  for  VWN,  3.42  eV  for 
HL)  is  considerably  larger  than  the  experimental  value10  (2.7  eV).  It  appears  to 
be  a  common  failing  of  local  density  functional  theory  that  it  gives  too  much 
binding1,  5-  20  the  reason  being  that  the  large  underestimate  of  the  total  energies 
of  the  free  atom  and  of  the  solid  (4.3%  for  H,  0.1%  for  Ni)  due  to  the  local  den¬ 
sity  approximation  do  not  quite  cancel,  the  error  being  slightly  larger  for  the 
atom  than  for  the  (more  homogeneous)  solid. 

The  bridge  site  minimum  lies  only  0.08  eV  higher  than  the  center  site 
minimum  and  the  top  site  minimum  lies  about  0.32  eV  higher.  It  ought  to  be 
mentioned  here  that  although  the  individual  curves  have  a  relative  numerical  pre¬ 
cision  of  better  than  0.01  eV,  there  may  be  a  relative  error  between  the  energies 
of  the  center  and  bridge  site  curves  of  0.1  eV  and  of  the  center  and  top  site 
curves  of  0.2  eV.  The  reason  for  this  is  that  since  H  atoms  approach  the  surface 
Ni  atoms  more  closely  at  the  bridge  and  top  sites,  one  is  forced  to  choose  smaller 
muffin-tin  radii  when  calculating  these  curves,  resulting  in  a  less  well  converged 
basis  set  and  greater  uncertainties  in  the  extrapolation  to  the  infinite  basis  size 
energies.  The  extrapolation  procedure  is  explained  in  greater  detail  later. 
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The  near  degeneracy  of  the  center  and  bridge  site  minima  implies  that  the  H 
atoms  will  have  large  rms  displacements  from  their  equilibrium  positions  or  may 
even  be  mobile  on  the  surface.  There  is  some  experimental  evidence  that  the  bar¬ 
rier  to  surface  diffusion  of  the  H  is  small.  Ion-channeling  experiments12  show 
that  the  H  atoms  have  a  rms  displacement  of  0.4  a,,  from  their  equilibrium  posi¬ 
tion.  Laser-induced  desorption  experiments21  show  that  the  barrier  to  diffusion  is 
about  0.17  eV  for  H  on  a  laser  treated  Ni(100)  surface.  Also  H  has  been  observed 
to  be  mobile  on  other  transition  metal  surfaces  by  field-ion  microscopy.22 

Also  shown  in  Fig.  1  are  other  theoretical  results  obtained  using  the  general¬ 
ized  valence  bond  (GVB)23  and  the  effective  medium24  methods.  The  inclusion  of 
configuration  effects  in  GVB  restricts  its  application  to  small  clusters  of  atoms 
which  may  not  adequately  represent  an  extended  surface  and  may  not  represent 
the  neighborhood  of  the  center,  bridge  and  top  sites  equally  well.  Also,  the  d 
electrons  are  not  included  explicitly.  The  effective  medium  technique  is  a  simple 
method  based  on  density  functional  theory.  It  has  the  advantage  that  it  is  fast 
enough  to  use  in  complicated  geometries  where  the  more  accurate  methods  used 
in  our  work  are  not  feasible  and  the  disadvantage  that  it  is  approximate,  in  par¬ 
ticular  the  interaction  of  the  adsorbate  with  the  metal  d-electrons  is  treated  in  a 


crude  way.  The  GVB  calculation  gives  an  equilibrium  position  in  good  agreement 
with  ours  but  the  relative  energy  difference  of  the  minima  of  the  curves  for  the 
three  sites  is  much  larger  than  in  our  calculation.  The  effective  medium 
approach  gives  basically  the  same  energy  differences  between  the  minima  of  the 
three  curves  as  our  calculation,  and  bond  lengths  that  are  all  somewhat  larger 
than  ours.  Also  the  detailed  shape  of  the  curves  is  somewhat  different  from  ours. 
The  semi-empirical  embedded  atom  method25  gives  energy  differences  in  good 
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Elcctron  energy  loss  (EELS)  experiments11  for  a  full  layer  of  H  on  Ni(100) 
yield  a  vibration  frequency  of  78  meV.  A  precise  calculation  of  this  frequency 
would  entail  a  calculation  of  the  total  energy  at  all  points  in  space  in  the  neigh¬ 
bourhood  of  the  minimum  energy  surface.  Then  the  difference  of  the  energies  of 
the  two  lowest  eigenstates  for  the  motion  of  the  H  atom  as  a  whole  would  give 
the  EELS  vibration  frequency.  This  has  been  done  within  the  effective  medium 
approach26  but  is  difficult  to  do  in  our  approach.  However  an  approximate  esti¬ 
mate  can  be  obtained  from  the  second  derivative  of  the  energy  curve  near  the 
minimum.  This  gives  a  value  of  90  meV  which  compares  reasonably  well  with 
experiment  given  the  approximations  we  have  made. 

Since  the  energies  of  interest  are  obtained  by  subtracting  total  energies 
which  are  several  orders  of  magnitude  larger  than  the  energies  of  interest,  it  is 
important  to  to  test  the  convergence  of  our  results  as  a  function  of  the  size  of  the 
various  expansion  sets  used,  and  to  test  the  sensitivity  of  our  results  to  the  choice 
of  various  parameters.  Here,  we  discuss  only  some  of  the  tests  performed,  namely 
convergence  with  respect  to  i)  the  size  of  the  basis  set  used  for  expanding  the 
wavefuncions,  and  ii)  the  number  of  layers  in  the  film,  and  sensitivity  to  the 
choice  of  the  exchange-correlation  potential. 

One  virtue  of  the  LAPW  method  is  that  the  basis  can  be  systematically 
improved  by  including  all  plane  waves  with  wavevectors  lying  within  successively 
larger  spheres  of  radius  K max.  Fig.  2a  shows  the  total  energy  versus  the  distance 
of  the  center  site  H  atoms  above  surface  Ni  layer  for  two  values  of 
where  RMT  is  the  muffin  tin  radius  of  the  center  layer  Ni  atoms.  In  Fig.  2b  the 
total  energy  for  r=0.6«9  is  plotted  versus  cip(-RMT  I\maJy  The  six  points 
shown  fit  the  exponential  with  an  rms  error  of  0.008  eV.  An  almost  equally  good 
fit  can  be  obtained  to  a  power  law  (rms  error  0.011  eV)  with  a  exponent  « -12.5. 
The  estimate  of  Etot  for  a  fully  converged  basis  (plotted  as  the  bottom  curve  in 
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Fig.  2a)  is  obtained  assuming  the  exponential  dependence.  Assuming  power  law 
dependence  would  give  a  curve  uniformly  shifted  down  by  0.08  eV.  The  position 
of  the  minimum  and  the  shape  of  the  curves  for  RMT  Kmax  =  7.5  and  8.0  are 
almost  identical.  Hence  it  is  possible  to  obtain  the  correct  vertical  position 
without  having  a  fully  converged  basis  as  has  been  noticed  by  other  authors/ 
However,  since  the  center,  bridge,  and  top  site  curves  are  calculated  with 
different  values  of  and  are  therefore  not  equally  well  converged,  it  is  very 
important  to  perform  the  extrapolation  to  obtain  the  relative  energy  differences 
between  them. 

The  binding  energies  of  H  at  two  positions  above  the  center  site  were  recal¬ 
culated  using  a  thinner  slab  consisting  of  3  Ni  layers  with  a  layer  of  H  on  each 
side  (shown  as  square  points  in  Fig.  1).  The  close  agreement  between  these  ener¬ 
gies  and  those  for  the  fatter  slab  give  us  confidence  that  the  energies  are  con¬ 
verged  with  respect  to  the  thickness  of  the  slab.27 

Fig.  1  also  shows  that  using  the  HL  exchange-correlation  potential  instead  of 
the  VWN  potential  results  in  a  binding  energy  curve  almost  parallel  to  that  for 
the  VWN  potential  over  the  range  of  distances  shown  but  shifted  up  by  0.19  eV. 

We  have  shown  that  accurate  total  energy  calculations  can  be  used  to  deter¬ 
mine  the  geometry  of  adsorbates  on  surfaces. 

It  is  a  pleasure  to  acknowledge  many  useful  discussions  with  Anders  Carls- 
son,  Gayanath  Fernando,  Matti  Manninen,  Jens  Norskov,  Mats  Persson,  and 
Mark  Stiles  and  assistance  with  computing  by  the  Center  for  Materials  Science  of 
Los  Alamos  National  Laboratory.  This  research  was  funded  by  the  Office  of 
Naval  Research. 
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Figure  Captions 


Figure  1.  Binding  energy  calculated  by  various  authors  as  a  function  of  the 

vertical  distance  of  the  H  atoms  above  the  surface  Ni  layer,  in  the  center  ( _ ) 

bridge  ( _ )  and  atop  (...)  sites.  The  calculated  points  are  generally  indicated 

by  dots.  The  two  squares  show  that  the  results  for  WVN  5-layer  are  similar  to 
VWN  7- layer. 

Figure  2.  Convergence  of  the  total  energy  Et*  with  respect  to  the  basis 
size.  Fig.  (2a)  shows  curves  of  Ef  ot  as  a  function  of  the  vertical  distance  of  the  H 
atoms  above  the  surface  Ni  layer  for  A'max  equal  to  7.5  and  8.  is  the 

muffin-tin  radius  of  the  central  layer  Ni  atoms  and  /vmax  is  the  largest  wavevec- 
tor  included  in  the  basis  functions.  The  estimated  curve  for  an  infinite  basis  is 
also  shown.  Fig.(2b)  shows  the  procedure  used  to  estimate  Etot  for  an  infinite 
basis,  a  fit  =1.78 


